(N 

o 

(N 
C 

oo 



Amino-acid-dependent main-chain torsion-energy terms for protein systems 

Yoshitake Sakae 1,2 and Yuko Okamoto 2,3,4 

^Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, 
Aichi 444-8585, Japan 

^Department of Physics, Graduate School of Science, Nagoya University, Nagoya, Aichi 464-8602, 
Japan 

'^Structural Biology Research Center, Graduate School of Science, Nagoya University, Nagoya, Aichi 464-8602, 
Japan 

4 ' Center for Computational Science, Graduate School of Engineering, Nagoya University, Nagoya, Aichi 464-8603, 
Japan 

(Dated: 8 March 2013) 

Many commonly used force fields for protein systems such as AMBER, CHARMM, GROMACS, OPLS, and 
ECEPP have amino-acid-independent force-field parameters of main-chain torsion-energy terms. Here, we 
propose a new type of amino-acid-dependent torsion-energy terms in the force fields. As an example, we 
applied this approach to AMBER 003 force field and determined new amino-acid-dependent parameters for 
ib and vj' angles for each amino acid by using our optimization method, which is one of the knowledge- 
based approach. In order to test the validity of the new force-field parameters, we then performed folding 
simulations of a-helical and /3-hairpin peptides, using the optimized force field. The results showed that 
the new force-field parameters gave structures more consistent with the experimental implications than the 
original AMBER ff03 force field. 
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I. INTRODUCTION 

Computer simulations of protein folding into native 
structures can be achieved when both of the following 
two requirements are met: (1) potential energy func- 
tions (or, force fields) for the protein systems are suf- 
ficiently accurate and (2) sufficiently powerful conforma- 
tional sampling methods are available. Professor Harold 
A. Scheraga has been one of the most important pioneers 
in studies of both of the above requirements^. By the 
developments of the generalized-ensemble algorithms (for 
reviews, see, e.g., Refs!^) and related methods, Require- 
ment (2) seems to be almost fulfilled. In this article, we 
therefore concentrate our attention on Requirement (1). 

There are several well-known all-atom (or united- 
atom) force fields, such as AMBERS, CHARMMikil, 
OIM.S 12 1 \ GROMOS±i±£, GROMACS±&i£, and 
ECEP P 18,19 . Generally, the force-field parameters are 
determined based on experimental results for small 
molecules and theoretical results using quantum chem- 
istry calculations of small peptides such as alanine 
dipeptide. 

In a force field, the potential energy is usually com- 
posed of the bond-stretching term, the bond-bending 
term, the torsion-energy term, and the nonbonded en- 
ergy term. In these energy terms, it is known that the 
torsion-energy term is the most problematic. For in- 
stance, the ff94£ and ff96^ versions of AMBER differ 
only in the main-chain torsion-energy parameters. Nev- 
ertheless, the sccondary-structure-forming tendencies of 
the two force fields are quite different 2 ^— . Therefore, 
many researchers have studied this main-chain torsion- 
energy terms and their force-field parameters. For in- 
stance, newer force-field parameters of the main-chain 
torsion-energy terms about <f> and tp angles have been 
developed, which are, e.g., AMBER 599861, AMBER 



ff032, CHARMM22/CMAP— and OPLS-AA/LH. The 
methods of the force-field refinement thus mainly concen- 
trate on the torsion-energy terms. These modifications of 
the torsion energy are usually based on quantum chem- 
istry calculations^ 6 - - — or NMR experimental result a 31 ' 32 . 

We have also proposed a new main-chain torsion- 
energy term, which is represented by a double Fourier 
series in two variables, the main-chain dihedral angles <fi 
and ?/ ' 33 ' 34 . This expression gives a natural representa- 
tion of the torsion energy in the Ramachandran spaced 
in the sense that any two-dimensional energy surface pe- 
riodic in both (j> and vb can be expanded by the double 
Fourier series. We can then easily control secondary- 
structurc-forming tendencies by modifying the main- 
chain torsion-energy surface. We have presented pre- 
liminary results for AMBER ff94 and AMBER ff96^^. 
Moreover, we have introduced several optimization meth- 
ods of force-field parameters 2 ! ' 25 ' 36 ' 37 . These methods 
are based on the minimization of some score functions 
by simulations in the force-field parameter space, where 
the score functions are derived from the protein coordi- 
nate data in the Protein Data Bank (PDB). One of the 
score functions consists of the sum of the square of the 
force acting on each atom in the proteins with the struc- 
tures from the PDB^r— . Other score functions are taken 
from the root-mean-square deviations between the orig- 
inal PDB structures and the corresponding minimized 
structure s 36 ' 37 . 

In this article, we propose a new type of the main-chain 
torsion-energy terms for protein systems, which can have 
amino-acid-dependent force-field parameters. As an ex- 
ample of this formulation, we applied this approach to the 
AMBER ff03 force field and determined new amino-acid- 
dependent main-chain torsion-energy parameters for ip 
(N-C a -C-N) and vb' (C^-C Q -C-N) by using our optimiza- 
tion method in RefsS 3 . - — . 
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In section 2 the details of the new main-chain torsion- Fourier coefficient. Namely, the energy term -Etoraion has 

energy terms are given. In section 3 the results of ap- 7n($) an d V n (&) as force- field parameters, 

plications of the method to AMBER ff03 force field and We can further write the torsion-energy term as 
those of folding simulations of two peptides are presented. 

Section 4 is devoted to conclusions. (mc) p(SC) /„>, 

-^torsion = -^torsion + -^torsion > 



II. METHODS 

A. Amino-acid-dependent force-field parameters 

The existing force fields for protein systems such as 
AMBERS, CHARMM±&±i, and OPLS 1 ^, etc. use 
essentially the same functional forms for the potential 
energy E con { except for minor differences. The confor- 
mational potential energy E con { can be written as, for 
instance, 

Econl = EbL + EbA + ^torsion + -Enonbond ■ (1) 

Here, E BL , E BA , E tmsion , and E nonhond represent 
the bond-stretching term, the bond-bending term, the 
torsion-energy term, and the nonbonded energy term, re- 
spectively. Each force field has similar but slightly dif- 
ferent parameter values. For example, the torsion-energy 
term is usually given by 

torsion = E E i 1 + C ° S ^ - 7n (*)]} , (2) 

where the first summation is taken over all dihedral an- 
gles $ (both in the main chain and in the side chains), 
n is the number of waves, j n is the phase, and V n is the 



where S tors idn an< ^ ^torsion are t ne torsion-energy terms 
for dihedral angles around main-chain bonds and around 
side-chain bonds, respectively. Examples of the dihe- 
dral angles in E^2n are <f> (C-N-C Q -C), V (N-C Q -C-N), 
0' (C^-Ca-N-C), V' (C r C Q -C-N), and w (C Q -C-N-C a ). 

(SC) 

The force-field parameters in i? torsion can readily depend 

on amino-acid residues. However, those in ^Msion are 
usually taken to be independent of amino-acid residues 
and the common parameter values are used for all the 
amino-acid residues (except for proline). This is because 
the amino-acid dependence of the force field is believed 
to be taken care of by the very existence of side chains. 
In Table I, we list examples of the parameter values for 
V> (N-C a -C-N) and ip' (C r C a -C-N) in general AMBER 
force fields. 

However, this amino-acid independence of the main- 
chain torsion-energy terms is not an absolute require- 
ment, because we are representing the entire force field 
by rather a small number of classical-mechanical terms. 
In order to reproduce the exact quantum-mechanical con- 
tributions, one can introduce amino-acid dependence on 
any force-field term including the main-chain torsion- 
energy terms. Hence, we can generalize E^^ n in Eq. (|3|) 
from the expression in Eq. @ to the following amino- 
acid-dependent form: 
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where k (= 1, 2, • ■ ■ , 20) is the label for the 20 kinds of 

amino-acid residues and are dihedral angles around 
the main-chain bonds in the fc-th amino-acid residue. 



B. Optimization method for force-field parameters 



amino-acid sequence homology as much as possible. If 
the force-field parameters are of ideal values, then all the 
chosen native structures are stable without any force act- 
ing on each atom in the molecules on the average. Hence, 
we expect 



In the previous subsection, we have generalized the where 
main-chain torsion-energy term E^^ on so that its pa- 
rameters are amino-acid dependent. The question is then 
how to obtain optimal parameter values for this new 
main-chain torsion-energy term. 

One method is to use the parameter optimization an d 
method that was introduced in Refs<22r— . We first re- 
trieve N native structures (one structure per protein) 
from PDB. We try to choose proteins from different 
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Here, N m is the total number of atoms in molecule m, 
i?toTai is the total potential energy for molecule m, and 
fi is the force acting on atom i. In reality, 0, and 
because F > 0, we can optimize the force-field parame- 
ters by minimizing F with respect to these parameters in 
the main-chain torsion-energy term in Eq. In prac- 
tice, we perform a minimization simulation in the main- 
chain torsion-energy force-field parameter space for this 
minimization^ 3 - - — . 



III. RESULTS AND DISCUSSION 

A. An example of the amino-acid-dependent force-field 
parameter optimizations 

We present the results of our optimizations of the 
force-field parameters Vi($^q) for the main-chain an- 
gles = ip {k) (N-C Q -C-N) and (C^-CVC-N) in 
Eq. flU). We did this for the case of AMBER ff03 force 

field. We determined these Vi(^>^ G ) values for the 19 
amino-acid residues except for proline. 

At first, we chose 100 PDB files with resolution 2.0 
A or better, with sequence similarity of amino acid 30.0 
% or lower, and with less than 200 residues (the aver- 
age number of residues is 117.0) from PDB-REPRDB2& 
(see Table HI)) . We then refined these selected 100 struc- 
tures. Generally, data from X-ray experiments do not 
have coordinates for hydrogen atoms. Therefore, we have 
to add hydrogen coordinates. Many protein simulation 
software packages provide with routines that add hydro- 
gen atoms to the PDB coordinates. We used the AM- 
BER11 program packaged. We thus minimized the total 
potential energy E tota i = E con{ + E solv + E constI with re- 
spect to the coordinates for each proten conformation, 
where -E^onstr is the harmonic constraint energy term 

(-Bconstr = Ehoavy atom ~ ^o) 2 ), and E so{v is the 

solvation energy term. Here, K x is the force constant of 
the restriction and xq are the original coordinate vectors 
of heavy atoms in PDB. As one can see from Econstr, the 
coordinates of hydrogen atoms will be mainly adjusted, 
but unnatural heavy-atom coordinates will also be mod- 
ified. We performed this minimization for all the 100 
protein structures separately and obtained 100 refined 
structures by using K x = 100 (kcal/mol). As for the 
solvation energy term E so \ v , we used the GB/SA solvent 
included in the AMBER program package (igb = 5 and 
gbsa = 1)*ML. 

For these refined protein structures, we performed the 
optimization of force-field parameters of ip and ip' 
angles for AMBER ff03 force field by using the fucn- 
tion F in Eq. ^ as the total potential energy function 
(-E-totai = E con f + -Esolv) for the Monte Carlo simulations 
in the parameter space. Here, we used AMBER11— for 
the force calculations in Eq. (O . We have to optimize the 
38 (= 2 x 19) parameters simultaneously by the simula- 
tions in 38 parameters. However, here, for simplicity, we 



just optimized two parameters, V\(ip^) and V\{ip'^)^ 
for each amino-acid residue k separately, keeping the 
other Vi values as the original values. In order to ob- 
tain the optimal parameters, we performed Monte Carlo 
simulations of two parameters (V± of ip and ip') for the 
19 amino-acid residues except for proline. In Table HTT1 
the optimized parameters are listed. 



B. Test simulations with two peptides 

In order to check the force-field parameters obtained 
by our optimization method, we performed the folding 
simulations using two peptides, namely, C-pcptidc of ri- 
bonuclease A and the C-terminal fragment of the Bl do- 
main of streptococcal protein G, which is sometimes re- 
ferred to as G-pcptidoS. The C-peptide has 13 residues 
and its amino-acid sequence is Lys-Glu _ -Thr-Ala-Ala- 
Ala-Lys+-Phc-Glu-Arg + -Gln-His + -Mct. This peptide 
has been extensively studied by experiments and is 
known to form an a-helix structur o 43 ' 44 . Because the 
charges at peptide termini are known to affect helix 
stabilit y 43 ' 44 , the N and C termini of the peptide was 
blocked with acetyl and N-mcthyl groups, respectively 
The G-peptide has 16 residues and its amino-acid se- 
quence is Gly-Glu--Trp-Thr-Tyr-Asp--Asp--Ala-Thr- 
Lys^Thr-Phe-Thr-Val-Thr-Glu 1 . The termini were 
kept as the usual zwitter ionic states, following the ex- 
perimental condition s 42 ' 45 ' 46 . This peptide is known to 
form a /3-hairpin structure by experiment s 42 ' 45,46 . 

For the folding simulations, we used replica-exchange 
molecular dynamics (REMD)i£. REMD is one of the 
generalized-ensemble algorithms, and has high conforma- 
tional sampling efficiency by allowing configurations to 
heat up and cool down while maintaining proper Boltz- 
mann distributions. We used the AMBER11 program 
package 3 ^. The unit time step was set to 2.0 fs, and 
the bonds involving hydrogen atoms were constrained by 
SHAKE algorithm 4 ^. Each simulation was carried out 
for 30.0 ns (hence, it consisted of 15,000,000 MD steps) 
with 16 replicas by using Langcvin dynamics. The ex- 
change procedure for each replica were performed every 
3,000 MD steps. The temperature was distributed ex- 
ponentially: 650, 612, 577, 544, 512, 483, 455, 428, 404, 
380, 358, 338, 318, 300, 282, and 266 K. As for solvent 
effects, we used the GB/SA model in the AMBER pro- 
gram package (igb — 5 and gbsa = l) 40 ' 41 . The initial 
conformations for each peptide were fully extended ones 
for all the replicas. The REMD simulations were per- 
formed with different sets of randomly generated initial 
velocities for each replica. 

In Fig. [TJ a-helicity and /3-strandness of two peptides 
obtained from the REMD simulations are shown. Wc 
checked the secondary-structure formations by using the 
DSSP program 4 ^, which is based on the formations of 
the intra-main-chain hydrogen bonds. As is shown in 
Fig. [TJ for the original AMBER ff03 force field, the a- 
hclicity is clearly higher than the /3-strandness not only 
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in C-peptide but also in G-peptide. Namely, the original 
AMBER ff03 force field clearly favors a-helix and does 
not favor /3-structure. On the other hand, for the opti- 
mized force field, in the case of C-peptide, the a-helicity 
is higher than the /3-strandness, and in the case of G- 
peptide, the /3-strandness is higher than the a-helicity. 
We conclude that these results obtained from the opti- 
mized force field are in better agreement with the experi- 
mental results in comparison with the original force field. 
In Fig. 3i -helicity and 7r-helicity of two peptides ob- 
tained from the REMD simulations are shown. For 3io 
hclicity, there is no large difference for both force fields 
in C-peptide, and in the case of G-peptide, the value of 
the optimized force field slightly decreases in compari- 
son with the original force field. 7r-helicity has almost 
no value in the both cases of the original and optimized 
force fields in two peptides. 

In Fig. [3l a-helicity and /3-strandness as functions 
of temperature for the two peptides obtained from the 
REMD simulations arc shown. For a-hclicity, the values 
of both force fields decrease gradually from low temper- 
ature to high temperature in the case of C-pcptide. On 
the other hand, in the case of G-peptide, there are small 
peaks at around 300 K and 358 K for the original and 
optimized force fields, respectively. For /3-strandness, in 
the case of C-pcptide, it is almost zero for both force 
fields. In the case of G-peptide, for the optimized force 
field, there is clearly a peak around 300 K. In Fig.[U 3io- 
hclicity and 7r-hclicity of the two peptides as functions 
of temperature arc shown. For 3io-hclicity, in the case 
of both peptides, the values of the optimized force field 
are lower than the original force field as a whole except 
around low temperature in C-peptide. For 7r-helicity, it 
is almost zero for both force fields in the two peptides. 

In Fig.[5j the lowest-energy conformations of C-peptide 
obtained from the REMD simulations in the case of the 
original and the optimized force fields are shown. In the 
case of the original force field, all the conformations have 
helices. No. 3 has only 3io-helix, No. 13 has both a- 
helix and 3io-helix, and the rest of the conformations 
have only a-helix. In the case of the optimized force 
field, seven conformations (Nos. 2, 4, 5, 7, 8, 12, and 13) 
have helices. Nos. 2, 5, 7, 13 have only a-helix, Nos. 4, 
12 have only 3io-helix, and No. 8 has both a-helix and 
3io-helix. Additionally, there is one /3-bridge structure 
in No. 10. In Fig. [51 the lowest-energy conformations 
of C-pcptide are shown. In the case of the original force 
field, all the conformations except for No. 4 have helices. 
No. 6 has both a-helix and 3io-helix, Nos. 5, 7, 11 have 
only 3io-helix, and the rest have only a-helix. In the case 
of the optimized force field, Nos. 11 and 12 have a-helix, 
No. 8 has 3io-helix, No. 5 has /3-bridge, and Nos. 7, 
9, 10 have /3-strand. These results clearly show that the 
optimized force field favors helix structure much less than 
the original force field, and, additionally, in the case of 
G-peptide, slightly favors /3-structure. 



IV. CONCLUSIONS 

The main-chain torsion-energy terms are the most 
problematic terms in the force field for protein systems. 
We therefore concentrate our attention on these terms in 
order to obtain optimal protein force field. In this article, 
we proposed amino-acid-dependent main-chain torsion- 
energy terms in the force field for protein systems. This 
generalization gives more freedom to the force-field opti- 
mization problem. In principle, we can introduce amino- 
acid dependence on any force-field term. The present 
work introduced this dependence on even the main-chain 
torsion-energy terms, which previously had been treated 
independent of the amino-acid residue type. 

As an example of the present general formalism, we 
modified the AMBER ff03 force field so that the Vi 
parameters of the main-chain ijj and ip' angles may be 
amino-acid dependent except for proline (hence, 38 pa- 
rameters were optimized). Although preliminary because 
we did not optimize the 38 parameters simulataneously, 
our optimized parameters already gave structures more 
consistent with the experimental implications than the 
original AMBER force field in the folding simulations of 
two small peptides. 

We can easily apply the present formulations to 
other popular force fields such as AMBER ff99SB, 
CHARMM22 / CMAP, etc. This will be our future work. 
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TABLE I. Torsion-energy parameters (V„ and 7 n ) for the 
main-chain dihedral angles ip and ip' in Eq. ([2]) for the origi- 
nal AMBER fT94, ff96, ff99, ff99SB, and ff03 force fields. The 
values are common among the amino-acid residues for each 
force field. Only the parameters for non-zero V n are listed. 



force field 




ip (N-CU-C-N) 






V (C 8 -C Q -C-N) 
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V n /2 


In 


n 


V n /2 


In 


ff94 
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0.75 


TV 


2 


0.07 







2 


1.35 


TT 


4 


0.10 







4 


0.40 


TT 








ff96 


1 


0.85 





2 


0.07 







2 


0.30 


TT 


4 


0.10 





ff99 


1 


1.70 


7T 


2 


0.07 







2 


2.00 


TT 


4 


0.10 





ff99SB 


1 


0.45 


7T 


1 


0.20 







2 


1.58 


TT 


2 


0.20 







3 


0.55 


TT 


3 


0.40 





ff03 


1 


0.6839 


7T 


1 


0.7784 


7T 




2 


1.4537 


TT 


2 


0.0657 


TT 




3 


0.4615 


TT 


3 


0.0560 






TABLE II. 100 proteins used in the optimization of force-field 
parameters. 



fold 


PDB ID 


chain 


PDB ID 


chain 


PDB ID 


chain 


PDB ID 


chain 


all a 


1DLW 


A 


1N1J 


B 


1U84 


A 


1HBK 


A 




1TX4 


A 


1V54 


E 


1SK7 


A 


1TQG 


A 




1V74 


B 


1DVO 


A 


1HFE 


S 


1J0P 


A 




1Y02 


A71-114 


1IJY 


A 


1I2T 


A 


1G8E 


A 




1VKE 


C 


1FS1 


A109-149 


1D9C 


A 


1AIL 


A 




1Q5Z 


A 


1T8K 


A 


10R7 


C 


1NG6 


A 




1C75 


A 


2LIS 


A 


1NH2 


B 


1Q2H 


A 




1NKP 


A 














all p 


1XAK 


A 


1T2W 


A 


1GMU 


Cl-70 


1AYO 


A 




1PK6 


A 


10FS 


B 


1BEH 


A 


1J08 


A 




1UXZ 


A 


1UB4 


C 


1LGP 


A 


1CQY 


A 




1PM4 


A 


10U8 


A 


1V76 


A 


1R6J 


A 




10A8 


D 


1IFG 


A 










a/p 


1IO0 


A 


1U7P 


A 


1JKE 


C 


1MXI 


A 




1LY1 


A 


1NRZ 


A 


1IM5 


A 


1VC1 


A 




lOGD 


A 


1IIB 


A 


1PYO 


D 


1MUG 


A 




1H75 


A 


1K66 


A 


1COZ 


A 


1D40 


A 


a + P 


1VCC 


A 


1PP0 


B 


1PZ4 


A 


1TU1 


A 




1Q2Y 


A 


1M4J 


A 


1N9L 


A 


1LQV 


B 




1A3A 


A 


1K2E 


A 


1TT8 


A 


1HUF 


A 




1SXR 


A 


1CYO 


A 


1ID0 


A 


1UCD 


A 




1F46 


B 


1KPF 


A 


1BYR 


A 


1Y60 


D 




1SEI 


A 


1RL6 


A 


1WM3 


A 


1FTH 


A 




1APY 


B 


1N13 


E 


1LTS 


C 


1UGI 


A 




1MWP 


A 


1PCF 


A 


1IHR 


B 


1H6H 


A 



TABLE III. Optimized V\/2 parameters for the main-chain 
dihedral angles if) and if>' for the 19 amino-acid residues (ex- 
cept for proline) in Eq. (j4|. The rest of the parameters are 
taken to be the same as in the original ff03 force field (see 
Table IJ). The original amino-acid-independent values are also 
listed for reference. 





lf> (N-C a -C-N) 


if)' (Cfl-Ca-C-N) 


original ff03 


0.6839 


0.7784 


Ala 


0.122 


0.150 


Arg 


0.409 


0.200 


Asn 


-0.074 


-0.162 


Asp 


-0.137 


0.182 


Cys 


0.361 


0.089 


Gin 


0.144 


-0.024 


Glu 


0.180 


0.152 


Gly 


0.258 




His 


0.020 


0.237 


Ik- 


0.643 


0.194 


Leu 


0.382 


0.257 


Lys 


0.222 


0.042 


Met 


0.141 


0.346 


Phe 


-0.010 


0.553 


Ser 


-0.248 


0.475 


Thr 


0.512 


0.328 


Trp 


0.027 


0.477 


Tyr 


0.082 


0.652 


Val 


0.142 


0.590 
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FIG. 1. a-helicity (a-1) and /3-strandness (a-2) of C-peptide 
and a-helicity (b-1) and /3-strandness (b-2) of G-peptide as 
functions of the residue number at 300 K. These values were 
obtained from the REMD simulations. Normal and dotted 
curves stand for the optimized and the original AMBER ff03 
force fields, respectivery. 
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FIG. 2. 3io-helicity (a-l) and 7r-helicity (a-2) of C-peptide 
and 3io-helicity (b-1) and 7r-helicity (b-2) of G-peptide as 
functions of the residue number at 300 K. These values were 
obtained from the REMD simulations. Normal and dotted 
curves stand for the optimized and the original AMBER ff03 
force fields, respectivery. 
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FIG. 3. a-helicity (a-1) and /3-strandness (a-2) of C-peptide 
and a-helicity (b-1) and /3-strandness (b-2) of G-peptide as 
functions of temperature. These values were obtained from 
the REMD simulations. Normal and dotted curves stand for 
the optimized and the original AMBER ff03 force fields, re- 
spectivery. 
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FIG. 4. 3io-helicity (a-1) and 7r-helicity (a-2) of C-peptide and 
3io-helicity (b-1) and 7r-helicity (b-2) of G-peptide as func- 
tions of temperature. These values were obtained from the 
REMD simulations. Normal and dotted curves stand for the 
optimized and the original AMBER ff03 force fields, respec- 
tivery. 



10 




FIG. 5. Lowest-energy conformations of C-peptide obtained 
for each replica from the REMD simulations, (a) and (b) are 
the results of the original AMBER ff03 and the optimized 
force fields, respectively. The conformations are ordered in 
the increasing order of energy. The figures were created with 
DS Visualizer' . 
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FIG. 6. Lowest-energy conformations of G-peptide obtained 
for each replica from the REMD simulations, (a) and (b) are 
the results of the original AMBER ff03 and the optimized 
force fields, respectively. The conformations are ordered in 
the increasing order of energy. The figures were created with 
DS Visualizer' . 



